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MENETELMA MONIKANAVAI SEN MITTAUS S IGNAALIN KAYTTAMI SEKS I 
LAHDEMALLINNUKSESSA 

KEKSINMON ALA 

Esilla. oleva keksinto liittyy uuteen ja kehit- 
tyneeseen ' menetelmaan monikanavaisen mittaussignaalin 
kayttamiseksi lahdemallinnuksessa. Erityisesti esilla 
oleva keksinto koskee uutta tapaa muuntaa pyorteetonta 
ja lahteetonta vektorikenttaa mittaavalla monikanavai- 
sella mittalaitteella mitatut mittaussignaalit lahdemal- 
linnuksen kannalta optimaaliseen muotoon. 

KEKSINNON TAUSTA 

Tutkittavan kohteen sahkoista toimintaa voidaan 
tarkastella esimerkiksi mittaamalla kohteen sahkovirto- 
jen tuottamaa tnagneettikenttaa kohteen ulkopuolelle si- 
joitetuilla antureilla. Nain saadun lahde j akauman mal- 
lintaminen mitattujen signaalien perusteella on kuiten- 
kin hyvin hankalaa, koska jokaisen magneett ikenttaj a- 
kauman tuottajana voi olla monta erilaista lahde j a- 
kaumaa. Toisin sanoen lahde jakaumaa ei voida ratkaista 
yksikasitteisesti mitattujen signaalien perusteella, jo- 
ten ongelman ratkaisemiseksi joudutaan asettamaan eri- 
laisia rajoitusehtoja, kuten joku etukateistietoon pe- 
rustuva parametrinen malli virralle tai parametrisoima- 
ton normirajoitus . 

Jatkuvan virt a j akauman parametrisoimattomaan 
mallintamiseen kaytetaan tavallisesti miniminormiesti- 
maatteja, joissa monikanavai sella mittalaitteella mitat- 
tu signaali pyritaan selittamaan virtajakaumalla, jonka 
normi on mahdollisimman pieni . Normiksi valitaan taval- 
lisesti LI- tai L2-normi, joista edellinen on virtaele- 
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menttien pituuksien summa yli valitun tilavuuden ja jal- 
kimmainen on virtaelementtien pituuksien nelioiden summa 
yli valitun tilavuuden. Miniminormiestimaattien lasken- 
taa on kuvattu mm. julkaisuissa "Interpreting magnetic 
5 fields of the brain: minimum norm estimates", M.S. 
Hamalainen et al. Medical & Biological Engineering & 
Computing, Vol. 32, pp. 35-42, 1994 seka ^^Visualization 
of magnetoencephalographic data using minimum current 
estimates", Uutela K. et al, Neurolmage, Vol. 10, pp. 

10 173-180, 1999. 

Perinteisten miniminormiestimaattien luontaisia 
ongelmia ovat laskennan raskaus ja kohinaherkkyys . Esi- 
merkiksi L2-normin tapauksessa tarvitaan kaanteismatrii- 
si matriisista G, jonka alkio (i,j) sisaltaa i:nnen ja 

15 j :nnen mittausanturin kytkentakenttien sisatulon, joten 
nama sisatulot joudutaan laskemaan kaikille anturipa- 
reille. Kytkentakentta maaritellaan siten, etta anturin 
mittaama signaali on virtaj akauman projektio kyseisen 
anturin kytkentakentalle . Kohinaongelmat johtuvat siita, 

20 etta antureille laskettu matriisi G on hairioaltis, jo- 
ten sen kaanteismatriisin laskennassa tarvitaan kaytan- 
non tilanteissa regularisoint ia . 

Regularisointimenetelmat , esimerkiksi singulaa- 
riarvoha j otelman katkaisuregularisoint i , ovat yleensa 

25 epaintuitiivisia, ja usein valttamatta tapauskohtaisesti 
ratkaistavia. Vaaranlainen regularisointi saattaa johtaa 
taysin virheelliseen mallinnustulokseen. 

• Nain ollen ongelmana lahdemallinnuksessa nykyi- 
sellaan on edelleen laskennan raskaus ja hitaus, mahdol- 

30 liset kohinan aiheuttamat virheet seka ' laskennan tapaus- 
kohtaisuus regularisoinnista johtuen. Edelleen, kuten 
ylla on.mainittu, regularisointi voi aiheuttaa merkitta- 
vasti virhetta lopulliseen laskentatulokseen. 
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KEKSINNON TARKOITUS 

Keksinnon tarkoituksena on poistaa edella mai- 
nitut epakohdat tai ainakin merkittavasti lieventaa nil- 
ta. Erityisesti keksinnon tarkoituksena on tuoda esiin 
5 uudentyyppinen menetelma, jolla voidaan merkittavasti 
keventaa ja nopeuttaa jatkuvan virtajakauman mallinnuk- 
seen liittyvaa laskentaa seka vah'entaa kohinaongelmia. 

Esilla olevan keksinnon tunnusomaisten piirtei- 
den osalta viitataan patenttivaatimuksiin. 
10 . 

KEKSINNON KUVAUS 

Esilla oleva keksinto liittyy uudenlaiseen ta- 
paan maarittaa tutkittavan kohteen jatkuva virtajakauma 
kayttaen mitatuista signaaleista laskettuja signaaliava- 

15 ruuden kantavektorikomponentteja. Kyseiset komponentit 
on valittu siten, etta ne kuvaavat tarkasteltavan virta- 
jakauman mahdollisimman riippumattomia piirteita, mika 
tehostaa ja tarkentaa laskentaa huomattavasti . 

Keksinnon perusidea on se, etta koska tavan- 

20 omaisella anturistolla anturikenttien sisatulojen lasku 
on raskasta ja hankalaa, kannattaa kayttaa kayttaa eri- 
koista anturistoa, jolla kytkentakentat ovat ortogonaa- 
lisia ja mielellaan analyyttisesti laskettavissa. .Tama 
voidaan periaatteessa toteuttaa sopivalla fysikaalisella 

25 anturistolla. Koska sopiva fysikaalinen anturisto on 
kuitenkin usein hankala valmistaa on useimmissa tapauk- 
sissa edullisempaa kayttaa tavanomaisesta sensoristosta 
laskennallisesti muodostettuj a virtuaalisensore j a , eli 
mittaussignaalit muunnetaan toisiksi sopivalla muunnok- 

30 sella siten, etta ne vastaavat signaaleja jotka virtuaa- 
linen mittauslaite olisi mitannut. Samalla voidaan tar- 
vittaessa poistaa ulkoisiin hairioihin liittyvat signaa- 
lit. Tama muunnos on kuvattu esimerkiksi patenttihake- 
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muksessa FI20030392, joka liitetaan tahan viittauksella. 
Muunnoksen jalkeen . lahdemallinnus suoritetaan optimaali- 
sella tavalla kayttaen signaaliavaruuden kantavektori- 
komponentteja varsinaisten mittaussignaalien sijaan. 
5 Eras keksinnon oleellinen piirre on se, etta muxinnoksen 
jalkeen lahdemallia ei enaa tarvitse regularisoida. 

Esilla olevan keksinnon kohteena on siis mene- 
telma kohteen virtaj akauman maarittamiseksi mittaamalla 
magneettikenttia kohteen laheisyydessa monikanavaisella 
10 mittalaitteella. Edullisesti kutakin kanavaa vastaa ai- 
nakin yksi mittausanturi ja kohdetta approksimoidaan 
pallosymmetrisella johteella. Kohde voi olla esimerkiksi 
ihmisen paa. 

Keksinnon mukaisesti muunnetaan monikanavainen 
15 mittaussignaali vastaten kutakin mittausanturia ennalta 
maaratyn virtuaalianturiston signaaleiksi ja maaritetaan 
tutkittavan kohteen virtaj akauma syvyydella r ennalta 
valitussa ortonormaalissa f unktiokannassa virtuaaliantu- 
riston signaaleista. Talloin virtaj akauman estimointi on 
20 nopeaa ja robust ia. Edelleen virtuaalianturistoa .vastaa- 
van signaalijoukon aikaansaamiseksi lasketaan moni- 
kanavaisesta mittaussignaalista multipoolikehitelma . Mul- 
tipoolikehitelma voidaan laskea kahdella tavalla; otta- 
malla huomioon tutkittavan kohteen ulkopuoliset magneetti- 
25 kentat tai huomioimatta niita. 

Edullisesti valitaan ortonormaaliksi funktiokanr 
naksi muotoa 

Af') = J^^Ci,„f(r)Xlm(e,(p) oleva kanta, 

/»0 w=-/ 

30 

missa f(jr) on vapaasti valittava radiaalinen 
funktio ja Xj^X^,(p) on ns , vektoripalloharmoninen. Talloin 

voidaan sijoittaa funktiokanta virtaj akaumayhtaioon ja 
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ratkaistaan virtaj akauman kertoimet analyyttisesti yhta- 
losta 

5 missa Yi on kertalukuun I liittyva vakio ja i? on 

tarkasteltavan pallon sade. Edullisesti funktiolla f(r) 
saadetaan virtajakaumamallin syvyyspainotusta. 

Lisaksi keksinnon kohteena on mittauslaite koh- 
-teen virtaj akauman maarittamiseksi mittaamalla magneet- 

10 tikenttia kohteen laheisyydessa . Mittauslaitteeseen kuu- 
luu joukko pyorteetonta ja lahteetonta vektorikenttaa 
tnittaavia mittauskanavia, jolloin kutakin kanavaa vastaa 
ainakin yksi mittausanturi , ja kasittelyvalineet mitta- 
ussignaalin kasittelemiseksi . Edullisesti kohdetta ap- 

15 proksimoidaan pallosymmetrisella johteella. 

Keksinnon mukaisesti kasittelyvalineisiin kuu- 
luu muuntomoduuli monikanavaisen mittaussignaalin vasta- 
ten kutakin mittausanturia muuntamiseksi • ennalta maara- 
tyn virtuaalianturiston signaaleiksi , ja laskentavali- 

20 neet tutkittavan kohteen virtaj akauman maarittamiseksi 
tai laskemiseksi syvyydella r ennalta valitussa ortonor- 
maalissa f unktiokannassa virtuaalianturiston signaaleis- 
ta. Eraassa sovelluksessa laskentavalineet on jarjestetty 
laskemaan monikanavaisesta mittaussignaalista miiltipooli- 

25 kehitelma. 

Keksinnon ansiosta jatkuvan virtaj akauman mal- 
linnukseen liittyvaa laskentaa voidaan merkittavasti ke-- 
ventaa ja nopeuttaa. Lisaksi keksinnon ansiosta voidaan 
vahentaa kohinaongelmia . Edelleen keksinto yksinkertais- 

30 taa lahdemallin regularisoinnin tai poistaa sen tarpeen, 
ja siten myos merkittavasti vahentaa virhemahdollisuut- 
ta. 



KUVALUETTELO 



Seuraavassa keksintoa selostetaan yksityiskoh- 
taisten esimerkkien avulla viittaatnalla oheiseen piirus- 
tukseen, jossa 

Kuvio 1 esittaa erasta esilla olevan keksinrion 
mukaista mitt'auslaitetta; ja 

Kuvio 2 esittaa vuokaaviota eraasta esilla ole- 
van keksinnon mukaisesta menetelmasta . 

KEKSINNON YKSITYISKOHTAINBN SELOSTUS 

Kuviossa 1 esitetaan eras esilla olevan keksin- 
non mukainen mittauslaite, johon kuuluu joukko pyortee- 
tonta ja lahteetonta vektorikenttaa mittaavia mittauska- 
navia l\ 1^, -./l^/ jolloin kutakin kanavaa vastaa ainakin 
yksi mittausanturi 2\ 2\ ...,2' ja kasittelyvalineet 3 
mittaussignaalin kasittelemiseksi . Edullisesti kasitte- 
lyvalineet on toteutettu tietokoneella. Edelleen kasitte- 
lyvalineisiin kuuluu muuntomoduuli 4 tnonikanavaisen mit- 
taussignaalin vastaten kutakin mittausanturia muuntami- 
seksi ennalta maaratyn virtuaalianturiston signaaleiksi 
ja laskentavalineet 5 tutkittavan kohteen virtajakauman 
maarittamiseksi . 

Edelleen kuviossa 1 on esitetty tutkittava koh- 
de K, jonka laheisyyteen mittausanturit 2\ 2% ...,2^ on 
jarjestetty. Tutkittavan kohteen sisalla on esitetty 
virtasilmukka, joka kuvaa mitattavan magneettikentan 
lahdetta. Tutkittava kohde voi olla ihmisen paa ja vir- 
talahteita voi olla useita. 

Kuviossa 2 esitetaan esilla olevan keksinnon 
eraan edullisen sovelluksen paavaiheet . Ensin anturis- 
tolla mitataan monikanavainen mittaussignaali , vaihe 21. 
Sen jalkeeri signaali muunnetaan vastaamaan ns. virtuaa- 
lianturistolla mitattua signaalia, vaihe 22, jolloin ma- 
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tetnaattiset operaatiot saadaan yksinkertaisemmiksi . Lo- 
puksi virtuaalianturiston signaalista voidaan yksinker- 
taisesti laskea virtajakauma kohteessa, vaihe 23, eli 
kaytannossa kuvata virtasilmukoiden paikat ja . voimakkuu- 
5 det pallon tai paan sisalla. 

Seuraayassa kuvataan keksinnoh matemaattinen 
tausta ja peruste. Kun magneettikentat muunnetaan Lapla- 
cen yhtalon perusratkaisuun r'^^^^^Ylm(0,(p) liittyviksi ker- 
toimiksi - a^^ +ibi„, , missa i on imaginaariyksikko, voi- 

10 daan ne lausua virta j akauman J(r) avulla pallokoordinaa- 
teissa {r^O^cp) , jolloin ne ovat muotoa 

^/., =7/ J^'^,« ••^(O^V (1), 



15 missa integrointi suoritetaan yli tutkittavan 

tilavuuden, 7, on kertalukuiin I liittyva vakio ja Xlm(0,q>) 

on ns, vektoripalloharmoninen, 

Tama muoto kay ilmi esim. julkaisusta '*Multipo- 
le expansions of electromagnetic fields using Debye po- 
20 tentials", C. G. Gray, American Journal of Physics , ' Vol . 
46, pp. 169-179, 1978. Edella mainittu lauseke on kyt- 
kentakenttamuotoa, jossa multipoolikertoimen M^^ kytken-- 

takentta on muotoa 
25 4,(0 = r'A^/„(e,9>) (2). 

Toisaalta vektoripalloharmoniset muodostavat 
syvyydella r ortonormaalin kannan, joten kyseisella sy- 
vyydella virtajakauma voidaan esittaa kyseisessa funk- 
30 tiokannassa : • 
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missa f (r) on joku radiaalinen funktio. 

Kun tarkasteltavaksi tilavuudeksi valitaan pal- 
lotilavuus, niin si joittamalla edellinen lauseke yhta- 
5 loon (1) virtajakauman kertoimet voidaan ratkaista ana- 
lyyttisesti: 

c,„=y,M,„,[j^%V(ryr] ' (4), 

10 missa y, on kertalukuim I liittyva vakio ja J? on 

tarkasteltavan pallon sade. Edellinen yhtalo (4) osoit- 
taa, etta ortonormaalissa kannassa esitetyn virtaja- 
kaumamallin kertoimet saadaan ratkaistua kertoimista M, 
taysin triviaalilla tavalla analyyttisin lausekkein il- 

15 man minkaanlaista regularisointia . Tama on laskennalli- 
sesti erittain nopeaa ja numeerisesti stabiilia. Funktio 
f(r) on vapaasti valittavissa ja silla voidaan saadella 
virtajakaumamallin syvyyspainotusta . 

Todettakoon viela, etta kantaf unktioina pallo- 

20 maisen johteen tapauksessa voi kayttaa muutakin ortogo- 
naalista kantaa tai kantaa, j ossa . sisatulot ovat muutoin 
nopeasti laskettavissa . Tallaisen saa esimerkiksi rikko- 
malla ortogonaalisuuden tunnetulla tavalla hiukan tai 
kayttamalla kantaa, joka ei ole ortogonaali, mutta jonka 

25 sisatulot ovat ennalta laskettavissa. 

Keksintoa ei rajata pelkastaan edella esitetty- 
ja sovellusesimerkkeja koskevaksi, vaan monet muunnokset 
ovat mahdollisia pysyttaessa patenttivaatimusten maarit- 
teleman keksinnollisen ajatuksen puitteissa. 



PATENTTIVAATIMUKSET 

1. Menetelma kohteen virtaj akauman maarittami- 
seksi mittaamalla magneettikenttia kohteen laheisyydessa 
pyorteetonta ja lahteetonta vektorikenttaa mittaavalla 

5 monikanavai sella mittalaitteella, jolloin kutakin kana- 
vaa vastaa ainakin yksi mittausanturi , tunnettu 
siita, etta 

muunnetaan monikanavainen mittaussignaali vas- 
taten kutakin mittausanturia -ennalta maaratyn virtuaa- 
10 lianturiston signaaleiksi, ja 

maaritetaan tutkittavan kohteen virtaj akauma 
ennalta valitussa tehokkaasti laskettavassa funktiokan- 
nassa virtuaalianturiston signaaleista . 

2. Patenttivaatimuksen 1 mukainen menetelma, 
15 tunnettu siita, etta kohdetta approksimoidaan pal- 

losymmetrisella johteella ja lasketaan monikanavaisesta 
mittaussignaalista kentan multipoolikehitelma . 

3. Patenttivaatimuksen 2 mukainen menetelma, 
tunnettu siita, etta lasketaan multipoolikehitelma' 

20 ottamalla huomioon tutkittavan kohteen ulkopuoliset mag- 
nee ttikent at : 

4. Patenttivaatimuksen 2 mukainen menetelma, 
tunnettu siita, etta lasketaan multipoolikehitelma 
huomioimatta tutkittavan kohteen ulkopuolisia magneetti- 

25 kenttia. 

5. Patenttivaatimuksen 2 mukainen menetelma, 
tunnettu siita, etta ulkopuoliset hairiot poistetaan 
jollakin muulla tunnetulla hairioidenpoistomenetelmalla 
ennen muunnosta. 

30 6, Patenttivaatimuksen 2 mukainen menetelma, 

tunnettu siita, etta valitaan ortonormaaliksi funk- 
tiokannaksi virtaj akaumayhtalo muotoa 
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missa f (r) on vapaasti valittava radiaalinen 
funktio ja Xlm(0,(p) on ns • vektoripalloharmoninen. 

7. Patenttivaatimuksen 4 mukainen menetelma, 
tunnettu siita, etta 

5 sijoitetaan ortonormaali funktiokanta virtaja- 

kaumayhtaloon, ja . 

ratkaistaan virtajakauman kertoimet analyytti- 
sesti yhtalosta 

10 missa 7, on kertalukuun 1 liittyva vakio ja i? on 

tarkasteltavan pallon sade ja Xlm{e,q>) on ns. vektoripal- 
loharmoninen . 

8. Patenttivaatimuksen 4 mukainen menetelma, 
tunnettu siita, etta saadetaan funktiolla f (r) vir- 

15 tajakaumamallin syvyyspainotusta . 

9. Mittauslaite kohteen virtajakauman maaritta- 
miseksi mittaamalla magneett ikenttia kohteen laheisyy- 
dessa, johon mittauslaitteeseen kuuluu 

joukko pyorteetonta ja lahteetonta vektorikent- 
20 taa mittaavia mittauskanavia (1; l\ l\ l") , jolloin 
kutakin kanavaa vastaa ainakin yksi mittausanturi (2; 2\ 
2\ ...,2"), ja 

kasittelyvalineet (3) mittaussignaalin kasitte- 
lemiseksi ja jossa kohdetta approksimoidaan pallosymmet- 
25 risella johteella, tunnettu siita, etta 

kasittelyvalineisiin kuuluu muuntomoduuli (4) 
monikanavaisen mittaussignaalin vastaten kutakin mitta- 
usanturia muuntamiseksi ennalta maaratyn virtuaaliantu- 
riston signaaleiksi , ja 
30 , laskentavalineet (5) tutkittavan kohteen virta- 

jakauman maarittamiseksi syvyydella. r ennalta valitussa 
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ortonormaalissa f unktiokannassa virtuaalianturiston sig- 
naaleista. 

10. Patenttivaatimuksen 9 mukainen mittauslai- 
te, t u n n e t t u siita, etta laskentavalineet (5) on 

5 jarjestetty laskemaan monikanavaisesta mittaussignaalista 
multipoolikehitelma . 

11. Patenttivaatimuksen 10 mukainen mittauslai- 
te, tunnettu siita, etta lasketaan multipoolikehi- 
telma ottamalla huomioon tutkittavan kohteen ulkopuoliset 

10 magneettikentat . 

12. Patenttivaatimuksen 10 mukainen mittauslai- 
te, tunnettu siita, etta lasketaan multipoolikehi- 
telma huomioimatta tutkittavan kohteen ulkopuolisia mag- 
neettikenttia . 

15 13. Patenttivaatimuksen 10 mukainen mittauslai- 

te, tunnettu siita, etta valitaan ortonormaaliksi 
funktiokannaksi virtajakaumayhtalo muotoa 

missa f (r) on vapaasti -valittava radiaalinen 

20 funktio. 

14. Patenttivaatimuksen 12 mukainen mittauslai- 
te, tunnettu siita, etta 

sijoitetaah ortonormaali funktiokanta virtaja- 
kaumayhtalopn, ja 
25 ratkaistaan virtajakauman kertoimet analyytti- 

sesti yhtalSsta 

missa on kertalukuun I liittyva vakio ja i? on 
tarkasteltavan pallon sade. 
30 15. Patenttivaatimuksen 13 mukainen mittauslai- 

te, tunnettu siita, etta saadetaan funktiolla f(r) 
virtajakaumamallin syvyyspainotusta. 
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16. Vaatimuksen 9 mukainen mittalaite ja ana- 
lyysiohjelmisto, missa mittalaite muuntaa signaalit virtu- 
aalisensoristoon ennen niiden tallennusta, ja analyysioh- 
jelmisto muuntaa tallennetun datan virtajakaumaksi . 



TIIVISTELMA 

Keksinnon kohteena on menetelma tulkita 
tutkittavan kohteen jatkuva virta j akauma 
kayttaen mitatuista signaaleista laskettu- 
ja kantavektorikomponentteja • Kyseiset 
komponentit on valittu siten, etta ne ku- 
vaavat tarkasteltavan virta jakauman mah- 
dollisimman riippumattomia piirteita, mika 
tehostaa ja tarkentaa laskentaa huomatta- 
vasti, Tama aikaansaadaan Tnuuntamalla mi- 
tatut signaalit virta jakauman kannalta 
luonnollisempaan muotoon ja poistaen sa- 
malla ulkoisiin hairioihin liittyvat sig- 
naalit kokonaan. Tallainen muunnos on ku- 
vattu esimerkiksi patent tihakemuksessa 
FI20030392. Muunnoksen jalkeen lahdemal- 
linnus suoritetaan optimaalisella tavalla 
kayttaen signaaliavaruuden kantavektori- 
komponentteja varsinaisten mittaussignaa- 
lien sijaan. Eras keksinnon oleellinen 
piirre on se, etta muunnoksen jalkeen lah- 
demallia ei enaa tarvitse regularisoida . 



(Fig 1) 
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Fig. 1 



Mitataan monikana- 
vaineri signaali 



Muunnetaan signaali 
virtuaalianturistoon 



Maaritetaan virta- 
jakauma 
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Fig. 2 
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